## [1] "ATTRIBUTE"  "WETLAND_TY" "ACRES"      "Shape_Leng" "Shape_Area"

## # A tibble: 9 × 2
##   hpu   result
##   <fct>  <dbl>
## 1 Bh    0.152 
## 2 Bhs   0.0472
## 3 Bhsm  0.0518
## 4 Bi    0.0790
## 5 E     0.0584
## 6 H     0.473 
## 7 I     0.181 
## 8 O     0.0285
## 9 T     0.0636
## Warning: Removed 4 rows containing non-finite values (`stat_ydensity()`).
## Warning: Groups with fewer than two data points have been dropped.
## Warning: Removed 4 rows containing missing values (`geom_point()`).

Is there a model?

## 
## Call:
## lm(formula = WET ~ hpu, data = hb_pedons_ext)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.18141 -0.06360 -0.04722  0.00153  0.76040 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.15169    0.01923   7.890 4.01e-14 ***
## hpuBhs      -0.10447    0.03066  -3.407 0.000733 ***
## hpuBhsm     -0.09989    0.05472  -1.826 0.068774 .  
## hpuBi       -0.07272    0.02867  -2.536 0.011638 *  
## hpuE        -0.09325    0.03663  -2.546 0.011337 *  
## hpuH         0.32131    0.11615   2.766 0.005974 ** 
## hpuI         0.02972    0.04374   0.679 0.497287    
## hpuO        -0.12319    0.04101  -3.004 0.002859 ** 
## hpuT        -0.08809    0.02494  -3.532 0.000468 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.162 on 346 degrees of freedom
##   (4 observations deleted due to missingness)
## Multiple R-squared:  0.09768,    Adjusted R-squared:  0.07682 
## F-statistic: 4.682 on 8 and 346 DF,  p-value: 1.855e-05
## Analysis of Variance Table
## 
## Response: WET
##            Df Sum Sq  Mean Sq F value    Pr(>F)    
## hpu         8 0.9830 0.122875  4.6822 1.855e-05 ***
## Residuals 346 9.0801 0.026243                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Warning: Removed 1 rows containing non-finite values (`stat_ydensity()`).
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 30 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 30 rows containing missing values (`geom_point()`).

Is there a model in W3?

## 
## Call:
## lm(formula = WET ~ hpu, data = w3_pedons_ext)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.34523 -0.04827 -0.01462  0.05837  0.24221 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)         0.12041    0.04312   2.792 0.009020 ** 
## hpuBhs podzol      -0.08286    0.06098  -1.359 0.184356    
## hpuBimodal podzol  -0.10821    0.09642  -1.122 0.270682    
## hpuE podzol        -0.05423    0.06312  -0.859 0.397090    
## hpuLithic Histosol  0.24317    0.06587   3.692 0.000884 ***
## hpuTypical podzol  -0.01292    0.06953  -0.186 0.853802    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.122 on 30 degrees of freedom
##   (1 observation deleted due to missingness)
## Multiple R-squared:  0.5043, Adjusted R-squared:  0.4217 
## F-statistic: 6.104 on 5 and 30 DF,  p-value: 0.0005175
## Analysis of Variance Table
## 
## Response: WET
##           Df  Sum Sq  Mean Sq F value    Pr(>F)    
## hpu        5 0.45405 0.090809  6.1045 0.0005175 ***
## Residuals 30 0.44627 0.014876                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Call:
## lm(formula = WET ~ water_table, data = w3_pedons_ext)
## 
## Residuals:
##        3        6       11       17       19       28       35 
## -0.10035 -0.06849 -0.08845 -0.01801  0.23055 -0.02270  0.06745 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)
## (Intercept)  0.185056   0.139867   1.323    0.243
## water_table -0.001216   0.002372  -0.513    0.630
## 
## Residual standard error: 0.1274 on 5 degrees of freedom
##   (30 observations deleted due to missingness)
## Multiple R-squared:  0.04992,    Adjusted R-squared:  -0.1401 
## F-statistic: 0.2627 on 1 and 5 DF,  p-value: 0.6301
## Analysis of Variance Table
## 
## Response: WET
##             Df   Sum Sq   Mean Sq F value Pr(>F)
## water_table  1 0.004263 0.0042626  0.2627 0.6301
## Residuals    5 0.081129 0.0162259

Watershed 3 wells

## Warning: package 'stringr' was built under R version 4.3.1

## Warning: Removed 109200 rows containing missing values (`geom_point()`).
## Warning: Removed 85993 rows containing missing values (`geom_point()`).
## Warning: Removed 272921 rows containing missing values (`geom_point()`).
## Warning: Removed 47973 rows containing missing values (`geom_point()`).
## Warning: Removed 442491 rows containing missing values (`geom_point()`).
## Warning: Removed 68058 rows containing missing values (`geom_point()`).
## Warning: Removed 85993 rows containing missing values (`geom_point()`).
## Warning: Removed 272921 rows containing missing values (`geom_point()`).
## Warning: Removed 199902 rows containing missing values (`geom_point()`).

Other soil moisture data

##    plot          WET
## 1   IL1 0.0133766429
## 2   IL2 0.0022351611
## 3   IL3 0.0056414600
## 4   IH1 0.0044546306
## 5   IH2 0.1531509824
## 6   E01 0.0064688146
## 7   E02 0.0006058321
## 8   E08 0.0008501202
## 9   E09 0.0064689358
## 10  E10 0.1881551786

## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 7674 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 19943 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 32704 rows containing non-finite values (`stat_smooth()`).
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 15780 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 32704 rows containing missing values (`geom_point()`).
## Warning: Removed 19943 rows containing missing values (`geom_point()`).
## Warning: Removed 7674 rows containing missing values (`geom_point()`).
## Warning: Removed 15780 rows containing missing values (`geom_point()`).

Old Denitrification soil moisture & oxygen

##   ID       UPL          WET
## 1  1 0.9997921 0.0002078759
## 2  2 0.9943766 0.0056234149
## 3  3 0.9409537 0.0590462370

## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Warning: Removed 26 rows containing missing values (`geom_point()`).
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Warning: Removed 87 rows containing missing values (`geom_point()`).
## Warning: Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).
## Removed 1 rows containing missing values (`geom_point()`).

## Warning: Removed 210 rows containing missing values (`geom_point()`).
## Warning: Removed 70 rows containing missing values (`geom_point()`).